Superspreading, overdispersion and their implications in the SARS-CoV-2 (COVID-19) pandemic: a systematic review and meta-analysis of the literature.

BACKGROUND: A recurrent feature of infectious diseases is the observation that different individuals show different levels of secondary transmission. This inter-individual variation in transmission potential is often quantified by the dispersion parameter k. Low values of k indicate a high degree of variability and a greater probability of superspreading events. Understanding k for COVID-19 across contexts can assist policy makers prepare for future pandemics. METHODS: A literature search following a systematic approach was carried out in PubMed, Embase, Web of Science, Cochrane Library, medRxiv, bioRxiv and arXiv to identify publications containing epidemiological findings on superspreading in COVID-19. Study characteristics, epidemiological data, including estimates for k and R0, and public health recommendations were extracted from relevant records. RESULTS: The literature search yielded 28 peer-reviewed studies. The mean k estimates ranged from 0.04 to 2.97. Among the 28 studies, 93% reported mean k estimates lower than one, which is considered as marked heterogeneity in inter-individual transmission potential. Recommended control measures were specifically aimed at preventing superspreading events. The combination of forward and backward contact tracing, timely confirmation of cases, rapid case isolation, vaccination and preventive measures were suggested as important components to suppress superspreading. CONCLUSIONS: Superspreading events were a major feature in the pandemic of SARS-CoV-2. On the one hand, this made outbreaks potentially more explosive but on the other hand also more responsive to public health interventions. Going forward, understanding k is critical for tailoring public health measures to high-risk groups and settings where superspreading events occur.

Exploring the role of superspreading events in SARS-CoV-2 outbreaks.

The novel coronavirus SARS-CoV-2 emerged in 2019 and subsequently spread throughout the world, causing over 600 million cases and 6 million deaths as of September 7th, 2022. Superspreading events (SSEs), defined here as public or social events that result in multiple infections over a short time span, have contributed to SARS-CoV-2 spread. In this work, we compare the dynamics of SSE-dominated SARS-CoV-2 outbreaks, defined here as outbreaks with relatively higher SSE rates, to the dynamics of non-SSE-dominated SARS-CoV-2 outbreaks. To accomplish this, we derive a continuous-time Markov chain (CTMC) SARS-CoV-2 model from an ordinary differential equation (ODE) SARS-CoV-2 model and incorporate SSEs using an events-based framework. We simulate our model under multiple scenarios using Gillespie's direct algorithm. The first scenario excludes hospitalization and quarantine; the second scenario includes hospitalization, quarantine, premature hospital discharge, and quarantine violation; and the third scenario includes hospitalization and quarantine but excludes premature hospital discharge and quarantine violation. We also vary quarantine violation rates. Results indicate that, with either no control or imperfect control, SSE-dominated outbreaks are more variable but less severe than non-SSE-dominated outbreaks, though the most severe SSE-dominated outbreaks are more severe than the most severe non-SSE-dominated outbreaks. We measure severity by the time it takes for 50 active infections to be achieved; more severe outbreaks do so more quickly. SSE-dominated outbreaks are also more sensitive to control measures, with premature hospital discharge and quarantine violation substantially reducing control measure effectiveness.

Spatial risk for a superspreading environment: Insights from six urban facilities in six global cities across four continents.

Introduction: This study sets out to provide scientific evidence on the spatial risk for the formation of a superspreading environment. Methods: Focusing on six common types of urban facilities (bars, cinemas, gyms and fitness centers, places of worship, public libraries and shopping malls), it first tests whether visitors' mobility characteristics differ systematically for different types of facility and at different locations. The study collects detailed human mobility and other locational data in Chicago, Hong Kong, London, São Paulo, Seoul and Zurich. Then, considering facility agglomeration, visitors' profile and the density of the population, facilities are classified into four potential spatial risk (PSR) classes. Finally, a kernel density function is employed to derive the risk surface in each city based on the spatial risk class and nature of activities. Results: Results of the human mobility analysis reflect the geographical and cultural context of various facilities, transport characteristics and people's lifestyle across cities. Consistent across the six global cities, geographical agglomeration is a risk factor for bars. For other urban facilities, the lack of agglomeration is a risk factor. Based on the spatial risk maps, some high-risk areas of superspreading are identified and discussed in each city. Discussion: Integrating activity-travel patterns in risk models can help identify areas that attract highly mobile visitors and are conducive to superspreading. Based on the findings, this study proposes a place-based strategy of non-pharmaceutical interventions that balance the control of the pandemic and the daily life of the urban population.

Quantifying Mycobacterium tuberculosis transmission dynamics across global settings: a systematic analysis.

The degree to which individual heterogeneity in the production of secondary cases ("superspreading") affects tuberculosis (TB) transmission has not been systematically studied. We searched for population-based or surveillance studies in which whole genome sequencing was used to estimate TB transmission and the size distributions of putative TB transmission clusters were enumerated. We fit cluster size distribution data to a negative binomial branching process model to jointly infer the transmission parameters $R$ (the reproductive number) and dispersion parameter, $k$, which quantifies the propensity of superspreading in a population (generally, lower values of $k$ ($<1.0$) suggest increased heterogeneity). Of 4,796 citations identified in our initial search, nine studies met inclusion criteria ($n=5$ all TB; $n=4$ drug resistant TB) from eight global settings. Estimated $R$ values (range: 0.10, 0.73) were below 1.0, consistent with declining epidemics in the included settings; estimated $k$ values were well below 1.0 (range: 0.02, 0.48), indicating the presence of substantial individual-level heterogeneity in transmission across all settings. We estimated that a minority of cases (range 2-31%) drive the majority (80%) of ongoing transmission at the population level. Identifying sources of heterogeneity and accounting for them in TB control may have a considerable impact on mitigating TB transmission.

Chaos of COVID-19 Superspreading Events: An Analysis Via a Data-driven Approach

Superspreading has become a key mechanism of COVID-19 transmission which creates chaos. The classical approach of compartmental models may not sufficiently reflect the epidemiological situation amid superspreading events (SSEs). We perform a data-driven approach and recognise the deterministic chaos of confirmed cases. The first derivative (≈difference of total confirmed cases) and the second derivative (≈difference of the first derivative) are used upon SSEs to showcase the chaos. Varying solution trajectories, sensitivity and numerical unpredictability are the chaotic characteristics discussed here. © 2023 Indian Institute of Health Management Research.

Extreme value modelling of SARS-CoV-2 community transmission using discrete generalized Pareto distributions.

Superspreading has been suggested to be a major driver of overall transmission in the case of SARS-CoV-2. It is, therefore, important to statistically investigate the tail features of superspreading events (SSEs) to better understand virus propagation and control. Our extreme value analysis of different sources of secondary case data indicates that case numbers of SSEs associated with SARS-CoV-2 may be fat-tailed, although substantially less so than predicted recently in the literature, but also less important relative to SSEs associated with SARS-CoV. The results caution against pooling data from both coronaviruses. This could provide policy- and decision-makers with a more reliable assessment of the tail exposure to SARS-CoV-2 contamination. Going further, we consider the broader problem of large community transmission. We study the tail behaviour of SARS-CoV-2 cluster cases documented both in official reports and in the media. Our results suggest that the observed cluster sizes have been fat-tailed in the vast majority of surveyed countries. We also give estimates and confidence intervals of the extreme potential risk for those countries. A key component of our methodology is up-to-date discrete generalized Pareto models which allow for maximum likelihood-based inference of data with a high degree of discreteness.

Dining-Out Behavior as a Proxy for the Superspreading Potential of SARS-CoV-2 Infections: Modeling Analysis.

BACKGROUND: While many studies evaluated the reliability of digital mobility metrics as a proxy of SARS-CoV-2 transmission potential, none examined the relationship between dining-out behavior and the superspreading potential of COVID-19. OBJECTIVE: We employed the mobility proxy of dining out in eateries to examine this association in Hong Kong with COVID-19 outbreaks highly characterized by superspreading events. METHODS: We retrieved the illness onset date and contact-tracing history of all laboratory-confirmed cases of COVID-19 from February 16, 2020, to April 30, 2021. We estimated the time-varying reproduction number (Rt) and dispersion parameter (k), a measure of superspreading potential, and related them to the mobility proxy of dining out in eateries. We compared the relative contribution to the superspreading potential with other common proxies derived by Google LLC and Apple Inc. RESULTS: A total of 6391 clusters involving 8375 cases were used in the estimation. A high correlation between dining-out mobility and superspreading potential was observed. Compared to other mobility proxies derived by Google and Apple, the mobility of dining-out behavior explained the highest variability of k (ΔR-sq=9.7%, 95% credible interval: 5.7% to 13.2%) and Rt (ΔR-sq=15.7%, 95% credible interval: 13.6% to 17.7%). CONCLUSIONS: We demonstrated that there was a strong link between dining-out behaviors and the superspreading potential of COVID-19. The methodological innovation suggests a further development using digital mobility proxies of dining-out patterns to generate early warnings of superspreading events.

Super-spreaders of novel coronaviruses that cause SARS, MERS and COVID-19: a systematic review.

PURPOSE: Most index cases with novel coronavirus infections transmit disease to just one or two other individuals, but some individuals "super-spread"-they infect many secondary cases. Understanding common factors that super-spreaders may share could inform outbreak models, and be used to guide contact tracing during outbreaks. METHODS: We searched in MEDLINE, Scopus, and preprints to identify studies about people documented as transmitting pathogens that cause SARS, MERS, or COVID-19 to at least nine other people. We extracted data to describe them by age, sex, location, occupation, activities, symptom severity, any underlying conditions, disease outcome and undertook quality assessment for outbreaks published by June 2021. RESULTS: The most typical super-spreader was a male age 40+. Most SARS or MERS super-spreaders were very symptomatic, the super-spreading occurred in hospital settings and frequently the individual died. In contrast, COVID-19 super-spreaders often had very mild disease and most COVID-19 super-spreading happened in community settings. CONCLUSIONS: SARS and MERS super-spreaders were often symptomatic, middle- or older-age adults who had a high mortality rate. In contrast, COVID-19 super-spreaders tended to have mild disease and were any adult age. More outbreak reports should be published with anonymized but useful demographic information to improve understanding of super-spreading, super-spreaders, and the settings in which super-spreading happens.

Erratum: Detection of SARS-CoV-2 infection clusters: The useful combination of spatiotemporal clustering and genomic analyses.

[This corrects the article DOI: 10.3389/fpubh.2022.1016169.].

Superspreading potentials of SARS-CoV-2 Delta variants across different contact settings in Eastern China: A retrospective observational study.

OBJECTIVES: As the genetic variants of SARS-CoV-2 continuously pose threats to global health, evaluating superspreading potentials of emerging genetic variants is of importance for region-wide control of COVID-19 outbreaks. METHODS: By using detailed epidemiological contact tracing data of test-positive COVID-19 cases collected between July and August 2021 in Nanjing and Yangzhou, China, we assessed the superspreading potential of outbreaks seeded by SARS-CoV-2 Delta variants. The transmission chains and case-clusters were constructed according to the individual-based surveillance data. We modelled the disease transmission as a classic branching process with transmission heterogeneity governed by negative binomial models. Subgroup analysis was conducted by different contact settings and age groups. RESULTS: We reported a considerable heterogeneity in the contact patterns and transmissibility of Delta variants in eastern China. We estimated an expected 14% (95% CI: 11-16%) of the most infectious cases generated 80% of the total transmission. CONCLUSIONS: Delta variants demonstrated a significant potential of superspreading under strict control measures and active COVID-19 detecting efforts. Enhancing the surveillance on disease transmissibility especially in high-risk settings, along with rapid contact tracing and case isolations would be one of the key factors to mitigate the epidemic caused by the emerging genetic variants of SARS-CoV-2.

Double impact: can our health care system deal with a terrorist attack during the COVID pandemic?

Health care organizations have been challenged by the COVID-19 pandemic since the first half of 2020. Both hospitals (especially emergency and intensive care departments) and ambulance services were overwhelmed by surging patient numbers during the 2 pandemic waves in 2020. In this study, the data of the 2016 multisite terrorist bombing attacks in Zaventem (Brussels International Airport) and Maalbeek (subway) are reviewed. It is simulated what the impact of similar attacks would be on an already challenged health care system and which COVID-19-specific measures would be favourable for the outcome. The limited access of ICU beds, operating rooms and surge capacity, as well as the number of COVID-positive victims are cardinal features challenging the medical response to mass casualty incidents of this magnitude. During the COVID-19 pandemic, disaster management is affected by the limited availability of intensive care beds and operation rooms, and the faltering reverse triage negatively influencing the response capacity. On the other hand, the impact of the COVID pandemic can also be favourable. Special concerns on a COVID-19-safe response are discussed. It must be avoided that the medical response and gathering of stranded passengers would become a superspreading event. Multisite terrorist attacks during a pandemic are possibly catastrophic for a health care system which is already beyond its limit in terms of surge capacity. COVID-19-specific recommendations for disaster management in case of terrorist attacks are provided.

Decrease in overdispersed secondary transmission of COVID-19 over time in Japan.

Coronavirus disease 2019 (COVID-19) has been described as having an overdispersed offspring distribution, i.e. high variation in the number of secondary transmissions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) per single primary COVID-19 case. Accordingly, countermeasures focused on high-risk settings and contact tracing could efficiently reduce secondary transmissions. However, as variants of concern with elevated transmissibility continue to emerge, controlling COVID-19 with such focused approaches has become difficult. It is vital to quantify temporal variations in the offspring distribution dispersibility. Here, we investigated offspring distributions for periods when the ancestral variant was still dominant (summer, 2020; wave 2) and when Alpha variant (B.1.1.7) was prevailing (spring, 2021; wave 4). The dispersion parameter (k) was estimated by analysing contact tracing data and fitting a negative binomial distribution to empirically observed offspring distributions from Nagano, Japan. The offspring distribution was less dispersed in wave 4 (k = 0.32; 95% confidence interval (CI) 0.24-0.43) than in wave 2 (k = 0.21 (95% CI 0.13-0.36)). A high proportion of household transmission was observed in wave 4, although the proportion of secondary transmissions generating more than five secondary cases did not vary over time. With this decreased variation, the effectiveness of risk group-focused interventions may be diminished.

A statistical framework for tracking the time-varying superspreading potential of COVID-19 epidemic.

Timely detection of an evolving event of an infectious disease with superspreading potential is imperative for territory-wide disease control as well as preventing future outbreaks. While the reproduction number (R) is a commonly-adopted metric for disease transmissibility, the transmission heterogeneity quantified by dispersion parameter k, a metric for superspreading potential is seldom tracked. In this study, we developed an estimation framework to track the time-varying risk of superspreading events (SSEs) and demonstrated the method using the three epidemic waves of COVID-19 in Hong Kong. Epidemiological contact tracing data of the confirmed COVID-19 cases from 23 January 2020 to 30 September 2021 were obtained. By applying branching process models, we jointly estimated the time-varying R and k. Individual-based outbreak simulations were conducted to compare the time-varying assessment of the superspreading potential with the typical non-time-varying estimate of k over a period of time. We found that the COVID-19 transmission in Hong Kong exhibited substantial superspreading during the initial phase of the epidemics, with only 1 % (95 % Credible interval [CrI]: 0.6-2 %), 5 % (95 % CrI: 3-7 %) and 10 % (95 % CrI: 8-14 %) of the most infectious cases generated 80 % of all transmission for the first, second and third epidemic waves, respectively. After implementing local public health interventions, R estimates dropped gradually and k estimates increased thereby reducing the risk of SSEs to approaching zero. Outbreak simulations indicated that the non-time-varying estimate of k may overlook the possibility of large outbreaks. Hence, an estimation of the time-varying k as a compliment of R as a monitoring of both disease transmissibility and superspreading potential, particularly when public health interventions were relaxed is crucial for minimizing the risk of future outbreaks.

Heterogeneity is a key factor describing the initial outbreak of COVID-19.

Assessing the transmission potential of emerging infectious diseases, such as COVID-19, is crucial for implementing prompt and effective intervention policies. The basic reproduction number is widely used to measure the severity of the early stages of disease outbreaks. The basic reproduction number of standard ordinary differential equation models is computed for homogeneous contact patterns; however, realistic contact patterns are far from homogeneous, specifically during the early stages of disease transmission. Heterogeneity of contact patterns can lead to superspreading events that show a significantly high level of heterogeneity in generating secondary infections. This is primarily due to the large variance in the contact patterns of complex human behaviours. Hence, in this work, we investigate the impacts of heterogeneity in contact patterns on the basic reproduction number by developing two distinct model frameworks: 1) an SEIR-Erlang ordinary differential equation model and 2) an SEIR stochastic agent-based model. Furthermore, we estimated the transmission probability of both models in the context of COVID-19 in South Korea. Our results highlighted the importance of heterogeneity in contact patterns and indicated that there should be more information than one quantity (the basic reproduction number as the mean quantity), such as a degree-specific basic reproduction number in the distributional sense when the contact pattern is highly heterogeneous.

Detection of SARS-CoV-2 infection clusters: The useful combination of spatiotemporal clustering and genomic analyses.

Background: The need for effective public health surveillance systems to track virus spread for targeted interventions was highlighted during the COVID-19 pandemic. It spurred an interest in the use of spatiotemporal clustering and genomic analyses to identify high-risk areas and track the spread of the SARS-CoV-2 virus. However, these two approaches are rarely combined in surveillance systems to complement each one's limitations; spatiotemporal clustering approaches usually consider only one source of virus transmission (i.e., the residential setting) to detect case clusters, while genomic studies require significant resources and processing time that can delay decision-making. Here, we clarify the differences and possible synergies of these two approaches in the context of infectious disease surveillance systems by investigating to what extent geographically-defined clusters are confirmed as transmission clusters based on genome sequences, and how genomic-based analyses can improve the epidemiological investigations associated with spatiotemporal cluster detection. Methods: For this purpose, we sequenced the SARS-CoV-2 genomes of 172 cases that were part of a collection of spatiotemporal clusters found in a Swiss state (Vaud) during the first epidemic wave. We subsequently examined intra-cluster genetic similarities and spatiotemporal distributions across virus genotypes. Results: Our results suggest that the congruence between the two approaches might depend on geographic features of the area (rural/urban) and epidemic context (e.g., lockdown). We also identified two potential superspreading events that started from cases in the main urban area of the state, leading to smaller spreading events in neighboring regions, as well as a large spreading in a geographically-isolated area. These superspreading events were characterized by specific mutations assumed to originate from Mulhouse and Milan, respectively. Our analyses propose synergistic benefits of using two complementary approaches in public health surveillance, saving resources and improving surveillance efficiency.

Differences in the superspreading potentials of COVID-19 across contact settings.

BACKGROUND: Superspreading events (SSEs) played a critical role in fueling the COVID-19 outbreaks. Although it is well-known that COVID-19 epidemics exhibited substantial superspreading potential, little is known about the risk of observing SSEs in different contact settings. In this study, we aimed to assess the potential of superspreading in different contact settings in Japan. METHOD: Transmission cluster data from Japan was collected between January and July 2020. Infector-infectee transmission pairs were constructed based on the contact tracing history. We fitted the data to negative binomial models to estimate the effective reproduction number (R) and dispersion parameter (k). Other epidemiological issues relating to the superspreading potential were also calculated. RESULTS: The overall estimated R and k are 0.561 (95% CrI: 0.496, 0.640) and 0.221 (95% CrI: 0.186, 0.262), respectively. The transmission in community, healthcare facilities and school manifest relatively higher superspreading potentials, compared to other contact settings. We inferred that 13.14% (95% CrI: 11.55%, 14.87%) of the most infectious cases generated 80% of the total transmission events. The probabilities of observing superspreading events for entire population and community, household, health care facilities, school, workplace contact settings are 1.75% (95% CrI: 1.57%, 1.99%), 0.49% (95% CrI: 0.22%, 1.18%), 0.07% (95% CrI: 0.06%, 0.08%), 0.67% (95% CrI: 0.31%, 1.21%), 0.33% (95% CrI: 0.13%, 0.94%), 0.32% (95% CrI: 0.21%, 0.60%), respectively. CONCLUSION: The different potentials of superspreading in contact settings highlighted the need to continuously monitoring the transmissibility accompanied with the dispersion parameter, to timely identify high risk settings favoring the occurrence of SSEs.

SARS-CoV-2 Transmission and Prevention in the Era of the Delta Variant.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) delta variant transmits much more rapidly than prior SARS-CoV-2 viruses. The primary mode of transmission is via short range aerosols that are emitted from the respiratory tract of an index case. There is marked heterogeneity in the spread of this virus, with 10% to 20% of index cases contributing to 80% of secondary cases, while most index cases have no subsequent transmissions. Vaccination, ventilation, masking, eye protection, and rapid case identification with contact tracing and isolation can all decrease the transmission of this virus.

Lockdowns exert selection pressure on overdispersion of SARS-CoV-2 variants.

The SARS-CoV-2 ancestral strain has caused pronounced superspreading events, reflecting a disease characterized by overdispersion, where about 10% of infected people cause 80% of infections. New variants of the disease have different person-to-person variability in viral load, suggesting for example that the Alpha (B.1.1.7) variant is more infectious but relatively less prone to superspreading. Meanwhile, non-pharmaceutical mitigation of the pandemic has focused on limiting social contacts (lockdowns, regulations on gatherings) and decreasing transmission risk through mask wearing and social distancing. Using a mathematical model, we show that the competitive advantage of disease variants may heavily depend on the restrictions imposed. In particular, we find that lockdowns exert an evolutionary pressure which favours variants with lower levels of overdispersion. Our results suggest that overdispersion is an evolutionarily unstable trait, with a tendency for more homogeneously spreading variants to eventually dominate.

Systematic review and meta-analyses of superspreading of SARS-CoV-2 infections.

Superspreading, or overdispersion in transmission, is a feature of SARS-CoV-2 transmission which results in surging epidemics and large clusters of infection. The dispersion parameter is a statistical parameter used to characterize and quantify heterogeneity. In the context of measuring transmissibility, it is analogous to measures of superspreading potential among populations by assuming that collective offspring distribution follows a negative-binomial distribution. We conducted a systematic review and meta-analysis on globally reported dispersion parameters of SARS-CoV-2 infection. All searches were carried out on 10 September 2021 in PubMed for articles published from 1 January 2020 to 10 September 2021. Multiple estimates of the dispersion parameter have been published for 17 studies, which could be related to where and when the data were obtained, in 8 countries (e.g. China, the United States, India, Indonesia, Israel, Japan, New Zealand and Singapore). High heterogeneity was reported among the included studies. The mean estimates of dispersion parameters range from 0.06 to 2.97 over eight countries, the pooled estimate was 0.55 (95% CI: 0.30, 0.79), with changing means over countries and decreasing slightly with the increasing reproduction number. The expected proportion of cases accounting for 80% of all transmissions is 19% (95% CrI: 7, 34) globally. The study location and method were found to be important drivers for diversity in estimates of dispersion parameters. While under high potential of superspreading, larger outbreaks could still occur with the import of the COVID-19 virus by traveling even when an epidemic seems to be under control.

Quantifying the importance and location of SARS-CoV-2 transmission events in large metropolitan areas.

Detailed characterization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission across different settings can help design less disruptive interventions. We used real-time, privacy-enhanced mobility data in the New York City, NY and Seattle, WA metropolitan areas to build a detailed agent-based model of SARS-CoV-2 infection to estimate the where, when, and magnitude of transmission events during the pandemic's first wave. We estimate that only 18% of individuals produce most infections (80%), with about 10% of events that can be considered superspreading events (SSEs). Although mass gatherings present an important risk for SSEs, we estimate that the bulk of transmission occurred in smaller events in settings like workplaces, grocery stores, or food venues. The places most important for transmission change during the pandemic and are different across cities, signaling the large underlying behavioral component underneath them. Our modeling complements case studies and epidemiological data and indicates that real-time tracking of transmission events could help evaluate and define targeted mitigation policies.